Singular matrix Darboux transformations in the 
inverse scattering method 



A. A. Pecheritsin, A. M. Pupasov and Boris F. Samsonov 

Physics Department, Tomsk State University, 634050 Tomsk, Russia. 
E-maih pech6r9phys.tsu.ru, pupasovaphys.tsu.ru, saiiisonov@phys.tsu.ru 



Abstract. Singular Darboux transformations, in contrast to the conventional 
ones, have a singular matrix as a coefficient before the derivative. We incorporated 
such transformations into a chain of conventional transformations and presented 
determinant formulas for the resulting action of the chain. A determinant 
representation of the Kohlhoff-von Geramb solution to the Marchenko equation 
is given. 



1. Introduction 

The idea to use Darboux (or equivalently supersymmetric, or simply SUSY) 
transformations for solving the inverse scattering problem for the one-dimensional 
one-channel Schrodinger equation was formulated in the most clear way by Sukumar 
[1]. It was essentially advanced by introducing phase equivalent transformations [2]- 
[5] generating potentials with modified spectra. This approach permitted to solve the 
long-standing problem of the deep or shallow nature of the nucleus-nucleus potentials 
[2] (see also the review paper [5]). In [7] this approach was reformulated by replacing a 
chain of phase equivalent (more precisely isophase) transformations by an equivalent 
A'^th-order transformation and using its Wronskian representation based on Crum- 
Krein formulas (HI IH] ■ In this way both the interaction potential and solutions of the 
corresponding Schrodinger equation acquire a compact analytic form which we call 
the determinant representation [7]. Moreover, in |10| the authors have shown how 
such SUSY transformations could yield correct phase shift effective range expansions. 

The generalization of this approach to multichannel scattering confronted with 
numerous difficulties. In particular, the problem of reconstructing, by the technique of 
Darboux transformations, the neutron-proton potential proposed in 1957 by Newton 
and Fulton was solved only quite recently [TT] . Obtained solution illustrates a general 
method of solving the multichannel inverse scattering problem which may become an 
alternative to the Gelfand-Levitan [12] or Marchenko [13] methods when applied to 
the two-channel inverse scattering problem. The main idea of the method is to divide 
the problem in two steps. This is possible since any 2x2 scattering matrix is uniquely 
determined by its eigenphase shifts and the mixing parameter which may be fitted to 
their experimental values independently. 

First, one constructs an uncoupled potential best fitting the given eigenphase 
shifts. Since the potential is supposed to be uncoupled, this problem reduces to 
solving n (n is the number of channels) "one-channel problems" . At this step, hence. 
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one can use well elaborated "one-channel techniques" (for a review see e.g. [B]). At 
the second step, using recently proposed eigenphase preserving (EPP) transformations 
|14| , one constructs the final potential by fitting the given mixing parameter only and 
keeping unchanged the eigenphase shifts. 

We have to emphasis that the formalism realizing one channel transformations 
is essentially different of that realizing two-channel EPP transformations. Of course, 
any first order in derivative one-channel transformation may be written in a matrix 
form but the rank of the matrix coefficient before the derivative is strictly less than n 
so that this matrix is non-invertible. We call such transformations singular Darboux 
(SUSY) transformations. By this reason the method developed in [15^ for replacing a 
chain containing such transformations by a single nth order transformation cannot be 
directly applied in this case. Note that such SUSY transformations, as far as we know, 
for the first time were studied in |16) . The main aim of the current paper is to extend 
the method of paper [TJ] by accepting singular transformations of a special type as 
links of Darboux transformations chain. This permits us to treat both singular and 
EPP transformations in a unified way. All intermediate quantities, such as solutions 
of the intermediate Schrodinger equations, do not enter into final expressions for the 
potential and wave functions so that they are expressed in terms of solutions of the 
initial equation only. 

The paper is organized as follows. In the next section we introduce necessary 
notations and describe the mathematical nature of the problem. In particular, 
we emphasis the difference between regular and singular Darboux transformations 
and present the latter in a form which allows us to treat both usual and singular 
transformations as links with equal rights of a transformation chain. The determinant 
formulas for the transformed solutions and potential are derived. In the third section 
we apply these results to construct the simplest 2x2 exactly-solvable potential 
describing '^S'l — ^Di neutron-proton elastic collisions. In Conclusion we discuss main 
results and outline a possible line of future investigations. 

2. Chains of Darboux transformations with singular links 

2.1. Preliminaries 

We start with the matrix Schrodinger equation 

Ho^E{x)^E'fE{x), Ho = -Ind^ + Vo{x), d := . (1) 

where e{x) = (?Ai£:(a;), . . . ipnE{x)y is a vector- valued function, /„ is n x n identity 
matrix, Vq{x) is annx n real and symmetric (potential) matrix. In general, x is a real 
variable which may belong to the whole real axis, to a semiaxis or to a closed interval 
of the real axis. In this section we do not associate equation ([1]) with any spectral 
problem and consider it as a system of differential equations of a special type. In the 
next section we shell treat it as a radial matrix Schrodinger equation. 

First, to fix notations, we would like to remind main notions about matrix 
Darboux transformations and their chains |15] . 

Suppose that we know N matrix solutions to equation ([T]) corresponding to 
different eigenvalue matrices A^ 7^ A; for fc 7^ 

HoUk=UkAk, k = l,...,N. 
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Here, in general, the spectral parameter A^. may be an arbitrary constant matrix 
but, following paper [13], we assume that it is a diagonal matrix with real entries. 
Moreover, for the case of equal thresholds, that we bear in mind, is a number. 

For the first transformation step we take matrix Ui (we call it the transformation 
matrix) and construct the transformation operator 

Li^o = Ind -F,, F, = idUi)U~' . (2) 

Note that it can be applied not only on vector- valued functions like but also on 
matrix- valued like ■ • • , i^N- In this way we get the matrix solutions V2 = Li^o^2, 
. . . , Vn — Li^qUn of the equation with the potential 

Vi=Va- 2dFi . 

Now V2 can be taken as transformation matrix for the Hamiltonian Hi — — + Vi 
to produce the potential 

V2 = Vi~ 2d{{dV2)V^^) =Vo- 2dF2 , = ^^1 + {dV2)V^^ 

and the transformation operator ^2^1 = Ind — {dV2)V2^ and so on, till one gets the 
potential 

Vn = Vo- 2dFM 

with Fn defined recursively 

FN^FN-i + {dYN)Yj;,\ N = l,2,... , Fo = 

and Yn being a matrix-valued solution to the Schrodinger equation at (N — l)th step 
of transformations, 

Hn-iYn — Yjsih-N , Hm-1 = —d^ + Vn-1 ■ 

The matrix Y^ is obtained by the action of the chain of — 1 transformations applied 
to the matrix Un , 

Yn = -^(^-1)^(^-2) • ■ • L24^iLi^qUn = i(Ar_i)^o^Af 

and it produces operator L^^(jv-i) = ^Ind + {dYiq)Yj^^ of the final transformation 
step for the chain of A^ transformations. 

To get in this way the final potential Vn resulting from the chain of A^ 
transformations, one has to calculate all intermediate transformation matrices Y}, 
j — 2, . . . , N performing a huge amount of unnecessary work even for the one-channel 
case. In practical calculations one is able to perform only few steps which restricts 
considerably possible applications of the method. Fortunately, for the one-channel 
case there exists what that are called Crum-Krein [51 [5] determinant formulas. Their 
multichannel generalization is given in [Tn|. These formulas allow one to omit all 
intermediate steps and go from Hq directly to H^. This is achieved by expressing the 
A^th order transformation operator 

Ln^q — Ln^{N-1) ■ ■ •-^2<-l-^1^0 

realizing the resulting action of the chain in terms of solutions of the initial Schrodinger 
equation ([T]) only. 

We would like to emphasis that the authors of [T3] considered chains of first order 
transformation operators of the form ^ where the coefficient before the derivative 
is a regular matrix which, being a;-independent, can always be reduced to identity 
matrix. Nevertheless, as it was first stressed by Andrianov et al [16], one of the main 



Singular matrix Darboux transformations in the inverse scattering method 



4 



features of matrix transformations is that the coefRcient before the derivative may be 
a singular matrix. The method developed in [T^ is not directly applicable to this case 
and needs some modifications. Below we give necessary modifications of that method 
for a particular type of singular matrices that appear when the SUSY method is used 
for solving the inverse scattering problem 



2.2. Singular matrix Darboux transformations of a special type 

Assume that the initial quantum system consists of two non-interacting subsystems 
(/) and (//) so that the matrix Vo{x) is block-diagonal 



Voix) 



Vo^'\x) 
V^"\x) 



Here Vq^^ is an m x m matrix, m=l,...,n— 1, and V^^^^"* is an {n — m) x [n — m) 
matrix. Since the total potential matrix Vo{x) is assumed to be real and symmetric, 
this implies that sub-matrices Vq'^^ and Vq^^^ are also real and symmetric. In this 
case, equation ([l} splits into two independent matrix equations 

^E'f^j^\x), hi'^ =-I„.d^ + V^'\x) (3) 

hi''H^^'\x) = E^f\x) , 4") ^ + , 

where *B^(a;) = (V'i£;(x), . . . , V'mfila;))*, ^''^^\x) = (?/>,m+i£; (a;), ...,■(/'«£; (a;))*. 

Assume that we want to realize the first order transformation over a subsystem 
only, say for definiteness over the subsystem (I). Evidently this should not affect the 
subsystem (II). Therefore the corresponding Darboux transformation operator, which 
we denote L^™\ has the form 

....... (^j- .1) w 

where the matrix U is an eigensolution of the Hamiltonian h\^^ 
h\^'^U = UK . 

We would like to emphasis that the coefficient before the derivative in ^ is 
a singular matrix. Therefore it can never be presented in the form Ind + w and 
such transformations cannot be directly incorporated into the usual chain of matrix 
Darboux transformations as considered in [15J . Nevertheless, theorems proven in [TS] 
have a more general character than Darboux transformation of the matrix Schrodinger 
equation. Actually, they represent a closure of a special recursion procedure. Below 
we show that with short modifications they may be applied to the current case also 
but first we will rewrite the operator Q in a form more suitable for our purpose. 

First we note that for any matrix U of the form 

u^(/i / )- (5) 

the following property holds 
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Next, if we move the matrix from the second term in the right hand side of 

equation (|3|) to its first term, the operator L^™) takes the form 

i^"^ ^ D,-a - {dU)U-\ 

Here the matrix vakied operator 

I-md 

Iji-^fj 



(6) 
(7) 



differentiates the first m entries of a vector <i> = {ipi, . . . , (pm, ^m+i, ■ ■ ■ ,Vn) and keeps 
unchanged any its component ipi for I > m. Similarly, while acting on an n x n 
matrix, it differentiates the first m entries of each column of this matrix and keeps 
unchanged any other entry. Below we will also use the composition of operators ([7]) , 
= Dm ■ ■ ■ ■ ■ Dm- When we need either differentiate or keep unchanged a separate 



component 1 <l < n oi vector $, we will use the notation 9(„i) defined as 

d(pu l<m, 
If I, I > m . 



a, 



(8) 



Note that the value of m is fixed by the dimension of the subsystem (/) . 

Expression ([6]) has the same structure as the usual Darboux transformation where 
Ind is replaced by Dm- It should also be noted that the matrix U is not solution to the 
initial Schrodinger equation ([T]) but it is constructed from solutions to the Schrodinger 
equation for the subsystem (/) ([3]). In the next section, we will present a generalization 
of the scheme developed in [15 which permits us to incorporate transformations such 
as the one given in ([6]) into a chain of usual transformations but first we need to 
introduce some new notations. 



2.3. Notations 

Let Uk, k = 1, . . . , N he a collection of n x n matrices 

\ 



Uk = 



U2,l;fc 



Ul,2;fe 
U2,2;fc 



U2' 



Any matrix Uk may be presented as either a collection of n-dimensional column- vectors 
Uj-k = {uij.k,---,Un,j-k)\ Uk = (J7i;fe, . . . , J7„;fe), A: = 1, . . . ,iV, or a collection of n- 
dimensional row-vectors Ul — • ■ • , Uj_n;k), j ~ ^, ■ ■ ■ ,n, Uk — {U^, . . . , UJ})*. 

Using these matrices, we define nM x nM matrix W 

( Ul U2 Us ... Um \ 

dUi DmU2 DmUs . . . DmUu 

d'Ui 92^2 DlMi ■■■ D-frMu ■ (9) 



W = W{Ui,U2,--.Mm) = 



1% 



D 



M-l 



Um ) 
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This matrix is used as the left upper block in a larger nN x uN, N > M, matrix 

Willi,...,UM;llM+i,..MN) 

( ly(M+l ■ ■ ■ Un \ 

W{UiM2,-.-Mm) DmUM+1 . . . DmUN 



'M 



D 



M-l 



dD^U 



N 



(10) 



Note that the first M matrices Uk, k = 1, . . . ,M in ([T(I| have the structure fixed by 
equation ([5]) whereas the structure of the matrices Uk for k = M + 1,. . . ,N is not 
fixed. For M — N all the matrices in this formula have the special structure ([5|) and 
the matrix (ITU)) reduces to 

Below we will also need a set of (nN + 1) x {nN + 1) matrices Wj, j = 1, . . . ,n 
obtained from (1101) by adding to the right a column constructed with the help of n- 
vector ^ = (V'l, . . . ,ip7i), to the bottom a row constructed from the jth row Ul of 
matrices Uk and the right bottom corner is filled with a derivative of jth element of 
the vector ^I*, 

WjiUi, . . . Mm'Mm+i,- . . ,Z^jv, *) = 
/ 



W{Ui 



,Um',Um+i, . . . ,Un) 



qN~M r,MTj3 



qN-MqM 



(11) 



(m) 



^.7 



Finally we introduce matrices Wi,j{Ui, . . . ,UmtUm+i, ■ ■ ■ ,Un), = l,...,n, 

obtained from the matrix (jlOp by replacing in its last matrix row the matrices 
QN'M-ij^M^^ (or D^-^Uk in the case N = M) by matrices Ul^ , k = 1, . . . , N , 

Wij{Ui, . . . ,Um',Um+i, ■ ■ ■ Mn) = 



W{UiM2.-.-Mm) 



DmUM+l 



d^'Ui d^'U2 

d^+^Ui d^'+^U2 



d 



N-2 



V 



Ul 



d^-^U2 



O^'Um 
O^'+'Um 



D^OUm+i 
OD^Um+i 



d^-'UM d 



N-AI-2 nM 



Dj^Um+i 



Un\ 

DmUN 

D^^^Un 
D^Un 
OD^Un 

:)N-M-2jjM^'^ 



(12) 



The matrices U]^-' 



i,j — l,...,n, k = 1,...,N, are constructed from the matrix 
gN-M-ijjM^^ (or D^-^Uk when M = A^) by replacing its jth row with the ith row 
of the matrix d^-'^^D^Ui (or dD^'^Uk when M = N). 
Before proving our main result we need two lemmas. 
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2-4-. Two lemmas 
Consider the matrix 

A = 



ap,i 
^1.1 

V ^2.1 



a 



0'l,p+n ^ 



b2,p+n / 



bl^p &l,p+l 
b2,p ^2,p+l 

Let a be the p x p submatrix of A with the entries Oij, i,j — 1, . . . ,p. Denote nijk 
the minor of A embordering a with jth {j — 1,2) row composed of bj^i, j = 1,2, 
i — 1, . . . ,p and fcth column {p < k < p + n). (For the definition of the embordering 
minor see Appendix.) Let also m*| be the minor obtained from rrijk by replacing 
its sth row composed of Osj, (s < p) with (p + t)th row composed of bt.i {t = 1,2). 
Let now a*" be obtained from a with the help of the same replacement, i.e. with the 
replacement of its sth row composed of a^.j, (s < p) by {p + i)th row of A composed 
of 6t,, (i = 1,2). 

Lemma 1. [15j // |a| 7^ then the following determinant identity takes place 



\a\m% 



The second lemma establishes a rule for differentiating the ratio of two 
determinants. 



Lemma 2. Let 



\WjiUi,...,UM;UM+l,--.MN,'^)\ 

\W{Ui,...Mm-Mm+u.-.Mn)\ 



(13) 



Then 



dipj = 



1=1 



Mn)\ 



\W{Ui,...,Um-Mm+u-.-Mn)\ 
Vj is the determinant of the matrix 
. ^U]\i;Um+i, ■ ■ ■ ,l^N, ^) by differentiating its last row. 



where 

Proof. In the expression for the derivative of the fraction ([T3l 

d\W{Ui,...MM-MM+l,...,UN)\ 



obtained from 



dip 



\W{Ui,...Mm-Mm+i,.--Mn)\ 
d\Wj(Ui,...MM-MM+i,.-.MN,^)\ 
\W{Ui,...,Um:.Um+i,---Mn)\ 



(14) 

we first analyze the derivative of the determinant \W{Ui, . . . ,IAm;Um+1t ■ ■ Mn)\- 
Taking into account the structure of matrix (jlO[) . one can see that a non zero 
contribution to this derivative gives a term corresponding to the differentiation of 
the last matrix row only 



d\W{Ui,...MM-MM 



■ Mn)V 



(15) 



1=1 



Note, that since Wi^i = for Z > m, then, in the case M 
expression goes up to / = m only. 



A^. the summation in this 
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Similarly, for the derivative of Wj , one gets 

n 

d\Wj{Ui, . . . ,Um;Um+i, ■ . ■ ,KN,'i')\ ^r, + (16) 

Here A^- is the determinant obtained from Wj by differentiating the Ith row in matrices 
of the next to last matrix row of Wj . 

We will use Lemma [T] to calculate A' . As matrix a we choose W. By adding 
to W two rows and one column, we obtain matrix A from Lemma [T] One row and 
the column are exactly the last row and the last column of the matrix Wj (fTTjl . The 
second added row is the ^th row of matrices d^^^'^'^'^D^~^Uk, k = 1,...,N, i.e. 
QN-M+ijjM~ijji^^ fc = 1, . . . , A^. Finally, to obtain matrix A we fill the right bottom 
corner with the element d^^'^^^^d^f^^ipi. Thus, identifying in Lemma [l]m*^ with 
A^-, a'^ with Wii, nijk with \Wj\, a^" with Wj,i and nitk with \Wi\, we find 

^ \WuiUi,...,UM;UM+i..---,UN)\\w,\ 

' \W[Uu...Mm-Mm+i,.-.Mn)\ 
\W,j{Ui, . . . Mm-Mm+i, ■ . ■ Mn)\\Wi\ 
\W{Ui,...Mm-Mm+i,...Mn)\ 
From here and it follows that 

d\w,{Uu ■ . . Mm-Mm+1,- --Mn, *)| = r, 

n 

WiA^i^ ■ ■ ■ ,Um;Um+i,- ■ ■ Mn)W]-\Wjj{Ui, . . . ,Um;Um+i,- ■ ■ ,Un)\vi] ■ 

We complete the proof by substituting (IT5|) and (fTTj) into (fT4|) . □ 

Note that this lemma plays a crucial role in the proof of our main results since 
it makes here applicable with only short modifications the proof of similar theorems 
given in [15]. 

2. 5. Transformation of a vector 

Let us consider a chain of N matrix first order Darboux transformations, where the 
first M transformations are singular (i.e., transformations of the subsystem (/)) and 
the remaining N — M transformations are usual matrix transformations. We will 
denote any particular transformation in this chain as j = \, . . . , N possibly 

with the superscript (m) if this is a singular transformation, L'^"^(^j_iy and by Lj^i^Q 
will be denoted the superposition of N first order transformations so that 

where 

L^^^k-i) ^ d - idYk)Y-\ k = M + l,...,N (18) 

and 

4"Vi) = - (^^fc)^"'' k = l,...M (19) 

with 



Fi^O, Yk^Lk^ol^k, fc = 2,...,7V. (20) 
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The resulting action of the chain is a transformation from the initial potential Vq to 
the final potential Vn- 

In the case M = N there are only singular transformations. 

We also note that the operator Ljv<-o can be applied not only on solutions of the 
Schrodinger equation but also on any vector- valued function ^/ = (-01, ... ,■(/)„ )*. 

Theorem 1. Chain of AI singular and N — M ordinary Darboux transformations acts 
on a vector ^/ ~ (-01, . . . ^ipnY o,s follows 

where 

\W,{Ui,...,Um-Mm+i,---Mn.-^)\ 



■ Mn)\ 



\W{Ui,...,Um-Mm+- 

with Wj{Ui,...MM'MM+i,---MNMN-,'^) and W{Ui,. 
defined in Ul\) and U(j\) respectively. 



Proof. We only give main ideas of the proof, since, at it was already mentioned, after 
proving Lemma [H the proof of the similar theorem given in is applicable here 
with short modifications. 

Following the method of paper 1151 , we will use the perfect induction method for 
proving the theorem. There are here two discrete variables, M and N. Therefore 
first we will prove the statement for M — N and then for a fixed value of M, we will 
comment on the case N > M. 

For M — N = 1 there is only one singular transformation 



Following the same lines as in [15^, one finds 



"Pi 



(21) 



(22) 



where 



\ dui. 



^1 \ 



meaning that the theorem holds for M = = 1. 

For M — N assume the theorem to hold for the chain of N 
transformations, i.e. 



^ _ \Wj{Uu.-.Mn-i,-^) 



\WiUi,...,UN-i)\ 
and prove it for the chain of N transformations. 
Operator L^^(jv-i) 



1 singular 

(23) 
(24) 



transforms the vector 8 into the vector $ 

L^^(-jv-i)6. To find the result of this transformation, i.e. the operator Lpf^Q 
ijv^(Ar-i)i(Ar_i)<_o (HZ]) it is ueccssary according to to calculate matrix Fjv 



L 



{N-iy 



-^qUn which is a collection of columns — (Yi- 



AT, ... , Yn-N). According to 
(|20p. each column transforms as an n- vector with the result given in (1241) where 5* 
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should be replaced by a column Ui-N of the matrix — {Ui-j^ . . 
for the entries yji of the vector Yj.pf = {yji, ■ ■ ■ l/jnY, one gets 

\W,{Uu...,UN-l,U^■,N)\ 



■Uu-nY- Therefore 



■ Mn- 



The determinant of the matrix Y^r — {yji ) can be calculated with the aid of the 
Sylvester (see Appendix) identity 

\W{Uu-.-Mn)\ 



\Y, 



\W{Uu...Mn^i)\ 



(25) 



According to (ITOl) , (|2ip and ([2^ , we can write down the components ipj of the vector 
$ as 



where 



yj 
^ N 



( 



Yf 



N 



(26) 



V • • • -yj" 

Note that in this case, the structure of matrix Yat coincides with the one given in ([5]), 
i.e., yji = Sji when i > m or j > m. Therefore the last row of matrix (1261) is zero 
for j > m except possibly for its right bottom corner and the first m elements when 
j < m. By this reason, it remains to calculate the non-zero entries of this row only 
and we will use Lemma [5] for that. 



Applying Lemma [2] to dyji yields 



dy. 



1 



\WiUi,...,l(N-l)\ 



1=1 



yH\W,,iiUi,...,UN-i)\ 



(27) 



Here Tji is the determinant of matrix Wj {Ui 



,Um'Mm 



+1) 



,Un-i, Ui-N) in which 



the elements of the last row are differentiated. Quite similarly, since for I < j < m 
the special operator ([5]) acts as the derivative, 9(TO)6'j — ddj, applying Lemma [2] to 
Pll) yields 



1 



\W{Ui,...,Un^i] 



Tj-Y,ei\Wj,i{Ui,...MN-i)\ 

1=1 



(28) 

..Mn-i,"^) by 



Here Tj is the determinant of the matrix obtained from WjiUi, 
differentiating its last row. 

From equations ((27|) and ((28|) it follows that for I < j < m the determinant 
of matrix Y^ equals the sum of two determinants. Furthermore, since the second 
terms of these formulas represent one and the same linear combination of yu in 
(|27p and 9i in (f28|) . which in their turn are elements of previous rows, the last 
row of the second determinant is a linear combination of the previous rows with 
the coefficients iM^jvl- ^^is reason this determinant vanishes. The matrix of the 
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first determinant consists of minors embordering the block WilAi, . . . ,Z-/jv-i) in matrix 
Wj{Ui, . . . Mni ^')- Applying the Sylvester identity to this determinant, we obtain 



When j > m, Dmdj = 9j and, dyji — because of the structure (O of matrix 

- N 



Therefore all elements of the last row of the matrix vanish except for the last one 



which is dj, and applying directly the Sylvester identity for calculating one gets 
the same result 



Thus, the determinant is given by ((29)) for all j = 1, . . . , n. 
Substituting this expression into (j26p and taking into account ((25)) . we get 

\W{U,,...Mn)\ ' ^^"^ 

thus finishing the proof of the theorem for M = N . 

For N ~ M, A/ + 1, . . ., the statement has just been proven for N ~ M. Moreover, 
any subsequent transformation is now the conventional Darboux transformation. 
Therefore, the proof in this case follows the same lines as in [l5\. □ 



2. 6. Transformation of potential 

Now we will show how to calculate the matrix potential of the Schrodinger equation 
obtained after TV Darboux transformations with M[< TV) singular transformations. 
Writing the transformed potential in the form 

Vn^Vq + 2dFN 

where 

FM^FN-i + {dYN)Yj,\ Fo^Q, Yn = L(^N~i)^oi^N , (31) 
we emphasis a recursive character of the procedure. 

Theorem 2. Let the matrix Fj^ be defined by the recursion i31\). Then the elements 
fi^j of the matrix _Fjv are expressed in terms of transformation matrices Uk, k ~ 
1, . . . , TV as follows 

^ \W{Ui,...Mm-Mm+i,...Mn)\ ■ 
where W{Ui, . . . ,Un) is defined by (ITOp and WijiJAi, . . . ,Un) is given in l\12^ . 

We do not dwell on the proof of the theorem since using Lemma[l]makes applicable 
with short modifications the existing proof of the similar theorem |15) . 



3. Determinant representation of the Kohlhoff-von Geramb solution of 
the Marchenko equation 

Note that the theorems from the previous section are not related to a spectral problem. 
They permit to construct matrix-valued potentials for which the multi-component 
Schrodinger equation, as a set of second order differential equations of a special type, 
can be solved exactly. In this section, as an illustration of this approach, we will apply 
obtained formulas for constructing the simplest local 2x2 potential matrix describing 
the ^Si — '^Di neutron-proton elastic scattering. Thus we will give a determinant 
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representation of the potential previously obtained by Kohlhoff and von Geramb [TS] 
who used Marchenko [T3] inversion method. 

Note that in this case x is the radial variable and we will use the conventional 
notation for it, x = r e (0, oo). 

We start with a diagonal matrix potential 

K,= (o g/^^), ^^€(0,00) 

describing the non-interacting free particle in s (subsystem (I)) and d (subsystem (II)) 
states. Thus we have here m = 1 and n ^ 2. 
Jost 

/.(fcr)=e^'=^ /,(fcr) = e^'^'- f 1 + - ^ 



and regular 

(fsikr) = i sin(fcr) , (fidikr) 



y kr [kr)"^ 
i [(3 — fc^r^) sin(fcr) — 3fcr cos(fcr)] 



one-channel s and d partial wave solutions are well known. 

We will now realize a four fold transformation over the composite system (I)+(II) 
choosing first two transformations singular with the transformation matrices hii, Ui 
having the structure (O 

U^^[^ <y5s(*fcir) ^ ^ ^ /s(-*fc2r) 

and two other transformations regular with the transformation matrices 

l>( -^s[K.r)N{-K) if,{Kr)N{K) \ _ , , 

Here iV(/c) = [{ki—ik){k2 — ik)]^^, k — xi^ + i), and x, fci, fc2 are free parameters. The 
factors N{±k) are introduced in (j32p to guaranty the symmetry, and hence, Hermitian 
character of the transformed potential. Evidently, if every intermediate potential 
emerging from any step of transformations is symmetric, the final potential is also 
symmetric. As it is shown in |17| . the matrix of the transformed potential remains 
to be symmetric after a first order transformation provided the self-Wronskian of the 
transformation matrix Yj defined as W[y, , Y,] = YjYj — (YJ)*Y, vanishes. The factors 
N{±k) are just chosen such that W[Y3, Y3] = W[F4, Y^] = 0. 

Note since we are not interested in intermediate potentials, the condition that they 
all are symmetric may be redundant. Actually, this is sufficient to impose a condition 
that the resulting action of the chain gives a symmetric potential. At present this 
problem is still waiting for its solution. 

Having fixed the set of matruces Uj, the potential matrix is constructed using 
Theorem 2, 

Vi{r) = Vir)^Voir)+AV{r), AV ^ (AV,,,) (33) 

with 

\ \w{u,,U2-MM)\ ) ^ ' 

where matrices Wij and W are given in ([TU| and ((T^ respectively. 
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Using the approach developed in |14) one can easily find the S'-matrix for the 
potential V ((M)) 

This is just the same S'-matrix that was used by Kohlhoff and von Geramb for 
constructing the simplest potential describing neutron-proton elastic collisions by the 
Marchenko inversion method. Thus we can state that the obtained potential V{r) 
([33]), (p4l) is a determinant representation of the Kohlhoff- von Geramb potential. 

We have to emphasis that according to property established in [M], EPP SUSY 
transformations change the order of channels. Therefore Vi,i and V2,2 correspond to 
d- and s-waves respectively. Also from ((33)) and (p4)) we extract central, Vc, tensor, 
Vt, and spin-orbital, Vq, components of the potential 

Vc = V2,2, Vt = Vi.2/V^, Vo^{V2a~Vxa/V2-Vi.i+&/r'')/i. (35) 

Choosing the same parameters fci , ^2 and x a-s in [TH [20] , we show these potential 
curves in Figure [TJ 

150 
100 
50 

> 

S 

-50 
-100 

12 3 

r, fm 

Figure 1. Exactly solvable potential (central Vc, tensor Vt and spin-orbit 
Vo terms) with parameters fci = 0.944, k2 = 0.232 fm^^ and x = 1-22 fm~^ 
taken from [20| . 

Comparing our results with those shown in figure 13 of paper |18) . one can see 
that in both cases the curves have the same shape. Note that the potentials ([55]) . ([M]) . 
as well as potentials constructed in |18) . are singular at the origin and by this reason 
they differ essentially from potentials obtained by Newton and Fulton [191 [20] . We 
believe that the difference with Ref. may be explained by a somewhat different 
input parameters. It is well-known that there is a family of potentials with the same 
scattering matrix, when there are bound states. In our case there is a single bound 
state and, hence, there could exist two additional parameters. Ad and As, known as 
the asymptotic normalization constants which provide an iso-phase deformation of the 
potential with the given scattering matrix. The values of these free parameters are 
fixed in both our model and the approach of Kohlhoff-von Geramb by the following 
reason. Physically acceptable components of a matrix potential should not have long- 
range tails. By this reason, corresponding asymptotic normalization constants can 
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be determined directly from the S'-matrix residue at the bound state pole |22] . In 
particular, the ratio of the asymptotic normalization constants reads 

_ Ad _ resS'2,i(fc = 1/02) 
^ As resS'i.i(fc — ik2) 

For the potential V{r) p3l) . ([M)) . as well as for the Kohlhofi-von Geramb potential, 
the parameter of 77 is fixed by the value 77 = fc|/(2x^) — 0.018081. 

4. Conclusion 

As it was pointed out by Andrianov et al [TB] , an essential feature of matrix Darboux 
transformations is that the coefficient before the derivative may be a singular matrix. 
We call such transformations singular transformations. Often these transformations 
allow one to find hidden symmetries of the problem |16| . Moreover, any one- 
channel transformation, being considered from a multichannel viewpoint, becomes 
singular. This leads to a problem of incorporating such transformations into a chain 
of conventional (i.e. regular) transformations. 

For the resulting action of the chain of conventional matrix Darboux 
transformations, there exists a generalization [T3] of the well known for the one-channel 
case Crum-Krein formulas [HIIS]- In the present paper, a method is developed which 
allows one to use one-channel and coupled-channel SUSY transformations on equal 
terms. We reformulated results obtained in [TS] such that singular transformations 
of a special type are included into the chain as links of the same rights as the 
conventional (regular) transformations except that they should be realized before the 
regular transformations. This approach together with recently introduced eigenphase 
preserving transformations |14j . may become an alternative to the Gelfand-Levitan- 
Marchenko method of solving the inverse scattering problem with the advantage that 
no needs exist for solving any integral equation. The values of the S'-matrix poles are 
incorporated into the the potential as parameters. 

We would like to emphasis that the use of exponential and spherical Bessel 
functions usually permits to fit the experimental scattering data with a very high 
precision. This means that both the potential and solutions of the corresponding 
Schrodinger equation are expressed via determinants (see theorems 1 and 2 above) 
containing elementary functions only. We believe that this is a big advantage as 
compared to the Gelfand-Levitan-Marchenko method. Moreover, in contrast to the 
usual scheme where one applies transformations step by step, our approach permits one 
to skip all intermediate calculations and obtain the final potential directly in terms 
of solutions of the initial Schrodinger equation, i.e. usually in terms of elementary 
functions. 

As an illustration of the method, using a special six-pole representation of the 
S matrix, we derived a determinant representation of the Kohlhoff-von Geramb [18] 
solution to the Marchenko equation. 

A drawback of the method, that we see, is that no a general recipe is known for 
choosing transformation matrices in a way to produce an Hermitian potential. Yet, in 
any concrete calculation this is enough to chose transformation matrices such that any 
intermediate potential is Hermitian. This is possible since corresponding condition is 
known [17j but the general problem is waiting for its solution. 
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Appendix A 

Here we formulate the Sylvester identity |21] . Consider a square matrix of dimension 
P + 9, P,? = 1,2, . . . 

^ ai_i . . . ai^p ai,p+i 
bi,i ■ ■ ■ bi^p 6i,p+i 

\ bq^l . . . bq^p bq^p+l 

Let a be the submatrix of dimension pxp composed of the elements atj, i, j = 1, . . . ,p. 
If to the bottom of a we add a line of elements bk,i, • . . , 6fe,p, to the right of a we add 
a column of elements ai,p+/, . . . , ap^p+i and the right bottom corner we fill with the 
element bk,p+i, we obtain a square matrix mjj. One says that rujj is obtained from A 
by embordering the block a with kth row and (p + Z)th column. The determinant 
is called an embordering minor in the determinant |y4|. Since k and I can take the 
values k,l = I, . . . ,q one has q x q embordering minors from which one can construct 
the matrix M = (mjj). The Sylvester identity relates the determinants \M\, \A\ and 
|a| as follows: 

\M\^\ar'\A\. (A.2) 
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